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We derive exact solitonic solutions of the one-dimensional time-dependent Gross-Pitaevskii equa- 
tion with time-dependent strengths of the harmonic external potential and the interatomic interac- 
tion. The time-dependence of the external potential and interatomic interaction are given in terms 
I of a general function of time. For an oscillating strength of the external potential, the solutions 

, correspond to breathing single and multiple solitons. The amplitude and frequency of the oscillat- 

ing potential can be used to control the dynamics of the center of mass of the solitons. For certain 
C ' values of these parameters, the solitons can be trapped at the center of the harmonic potential. 
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I. INTRODUCTION 

^ ' Since the experimental realization of solitons in Bose-Einstein condensates P, 0, S S S @| i intense interest in their 
properties has emerged 0, H, H, [13, HHj 113, [13 ■ Recently, exact solitonic solutions of the time-dependent Gross- 
Pitaevskii equation that describes the behavior of the condensate, have been obtained [H, [3, [I^, [H, [I3j [l3|- We 
have shown in that exact solutions of the Gross-Pitaevskii equation can not be obtained for general external 
potential and interatomic interaction. For the Gross-pitaevskii equation to be solved exactly, the external potential 
■ and the interatomic interaction must be related. As an example on such a relation, the external potential in the 
, Gross-Pitaevskii equation solved in Ref. [l3| is an expulsive harmonic potential and the interatomic interaction has an 
' exponential prefactor that grows with time with a rate that is related to the strength of the harmonic potential. Our 
O I investigation in 19] has provided a mathematical proof for the existence of such a correlation between the external 
, potential and the interatomic interaction. 

[ From an experimental point of view, it would be more interesting to find exact solutions for the external potentials 
T-H • and interatomic interactions used in the experiments. For instance, it would be more interesting to have exact 
^ ' solutions for the realistic case of an upright harmonic external potential and constant interatomic interaction rather 
, than for an inverted harmonic potential and exponentially growing interatomic interaction strength as in |16| . It 
turns out, however, that for the Gross-Pitaevskii equation to be solved exactly with an upright harmonic potential, 
the interatomic interaction will have to be oscillatory and divergent at regular times p9j . 

In this paper, we use the Darboux transformation method (2Qi] to investigate the possibility of obtaining exact 
^5 ', solutions of a time-dependent Gross-Pitaevskii equation with strengths of the harmonic external potential and inter- 
atomic interaction that are general functions of time. It turns out that some special cases of such a general form of 
harmonic potential and interatomic interaction can be easily realized experimentally by controlling the strength of 
the magnetic field that provides the harmonic trap and using Feshbach management of interatomic interaction [2]| . 
One interesting special case corresponds to an oscillating strength of the harmonic potential and practically constant 
strength of the interatomic interaction. The amplitude and frequency of the oscillating external potential can be used 
5J] to delay the escape of the soliton from the harmonic potential or even to trap it at the center. 

^ The rest of the paper is organized as follows. In the next section, we present the general form of the Gross-Pitaevskii 
equation to be solved. In section Hill we use the Darboux transformation method to derive the exact solutions. In 
section llVl we present and discuss the properties of the solutions in the special case of an oscillating external potential. 
We end in section |V] with a summary of our main conclusions. 



II. THE GROSS-PITAEVSKII EQUATION 



When the confinement of the Bose-Einstein condensate is much larger in, say the y and z (transverse) directions, 
compared to the confinement in the x direction, the system can be considered effectively one-dimensional along the 
X direction. The three-dimensional Gross-Pitaevskii equation can then be integrated over the y and z directions to 
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result in the following one-dimensional Gross-Pitaevskii equation [l9|, 

^dj^ = -^^1^ - \\p(t)x^i>{x,t) - 2ag(i)|V'(x,t)|V(x,t). (1) 

Here a is the s-wave scattering length, and A = IlOxI^i^, where uj^ and are the characteristic frequencies of the 
harmonic trapping potential in the x and transverse directions, respectively. In the last equation, length is scaled 
to aj^, time to and ip{x,t) to l/^/2a±, where a± = y^h/muj± is the characteristic length of the harmonic 

trap in the transverse direction. The dimensionless general functions p{t) and q{t) are introduced to account for the 
time-dependencies of the strengths of the trapping potential and the interatomic interaction. 

For Eq. [1] to be solved exactly using the Darboux transformation method, the functions p{t) and q{t) must be 
parametricUy related to each other throw a general function g{t) [Tlj. In this case, the Gross-Pitaevskii equation, 
that we derive exact solutions for, takes the form 

.aV^ = - [Xg{tr - m] - 2ae^-+^/^(')'^'|^(.,i)|V(x,t), (2) 

where ci is an arbitrary constant. 

III. DARBOUX TRANSFORMATION AND THE EXACT SOLUTIONS 

In the Darboux transformation method, we start by finding a linear system of equations for an auxiliary field ^'(a;, i) 
such that Eq. ^ is its consistency condition [H, Hoj. We find that Eq. ([2]) corresponds to the consistency condition 
of the following linear system [23j : 

*:,=CJ*A + P*, (3) 

*t = 2iC2j*A2 + 2<P*A + Aa;.gCJ*A + W*, (4) 

where, 

' ' \ (pi[x,t) (p2{x,t) y ' V ^1 
*^/^AiO\ ( ^Q[x,t) 

\0 \2 )' \ -^^Q*(x,^) 

_ f ia\Q{x,t)\'^ y/aXxg{t)Q{x,t) + iy/aQx{x,t) 

^ \ -^/aXxg{t)Q*{x,t) + iy/aQl{x,t) -ia\Q{x,t)\'^ 
and ({t) — expi^J Xg{t)dty 

The constants Ai and A2 are arbitrary constants. The subscripts x and t denote partial derivatives with respect to 
X and t, respectively. The function Q{x,t) is related to the wave function through Q{x,t) = ^{x,t)\/^e'''-~^^'^^^ 1^. 
Equation ^ is obtained from the consistency condition "^r^t — ^tx- 

The linear system of 8 equations, Eqs. ([3]) and reduces to an equivalent system of 4 equations with nontrivial 
solutions by making the following substitutions: Ai = — Aj, "02 =</>*, and ^2 — —^\- To be able to solve this reduced 
linear system we need to know an exact (seed) solution of Eq. Following our previous approach of finding seed 
solutions [III, we find the following seed solution to Eq. ^ 



'ipoix,t) = exp 



C2 -I- A.g/2 + i (^e^sfcoa: - ^° - ^Xgx^ + 2A^a J (e2(ci+c.+Ag) ^ ^2g2As^ 



(5) 



where C2, A and fco are arbitrary constants and g — dg/dt. 

The Darboux transformation can now be applied to the linear system to generate a new solution of Eq. ^ as 
follows 120] 

^ix,t)^Mx,t) + 2iX,+Xl)e-^^-''^^'/'+^'/'^^I^''j^^^ (6) 

Solving the linear system ^ and ^ using the seed solution Eq. ^ and then substituting for tl;o{x,t), ipi{x,t), and 
(f>i{x, t) in the last equation, we obtain the following new exact solution to Eq. 
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•1/4 —irix 

f] ' e ' 



{ 



-2i9i+ArVvx 



)(C393 



2iA^/ac4e 



I L2re2"''-'^'-v^^ + c2re-2"'?+A.v^:'= _ 4AciC2nA[(2Ai^ - A^)(cos26'i + cos6'4) 
+ (Ai - 2Air)(sin26'i + sin6'4) + (sin26'i + sin6'4)fco) |, 



(7) 



where 
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i\2A^a{2\T] - 1) + fco(fco(l - 2A?7) 
—2ar] — Ay/rjx, 



^/rjx) + 2Afco(Ai?7 + \/?7a;), 



Ai + fco + i(Ar + 2iAH - 2Air), 
A^Ai^ry - ((2Ah + h) ~ ^/ri)/^,/2, 
O2 = (fco(l - 2A?7) - 4AA^Air?7)/4A + 2A^a{\ - 2Ar/) + 2A,A(2Ai,ry - 
6*3 = (2v42a(2A?7 - 1) + ko{ko - 2Afco?7 + 2AV^a;))/2A, 
6*4 = (2AjAii - 2ArAir + Aifco)?? - V^A^x, 

= ko^/rjx + {2A^a - kl){2\r^ - 1)/2A, 
r = (A, - 2Ai,)2 + (A^ - 2Air)^ + ^A^a + fco(2A, - 4Ai, + fco), 
a = A^Aii + AiAir + Arfco/2, 
A^ = Re[v/(2Ai -ifco)2 -4A2a], 
A, = Im[V(2Ai-zfco)2-4A2a], 
7 = Ai + fco + i(Ar + 2iAii - 2Air), 

the subscripts r and i denote real and imaginary parts, respectively, and C3 and C4 are arbitrary constants. Further 
reduction and simplification of this new general solution as performed in [l3| will not be considered here. Furthermore, 
we do not attempt to obtain all classes of solitonic solutions. Instead, the focus in this paper will be on the effect of the 
parameters of the oscillating trapping potential on the dynamics of single and multiple solitons which are described 
by Eq. 0. 



IV. SPECIAL CASE: SOLITONS IN AN OSCILLATING HARMONIC TRAP 



In this section we derive from the general solution found in the previous section solitonic solutions for an oscillating 
harmonic trapping potential. To this end, we take 

r]{t) ^ j e"i+"2«'"M+'5)dt , (8) 

where ai, a2, w , and 5 are constants. Using the last equation and r/ — J e^^^dt to substitute for g in Eq. the 
Gross-Pitaevskii equation takes the form 

-2ae2'=i+("i+"^^'"'^*+*)/2|i/,(a;,i)p| (9) 

It should be noted that the specific choice of ■q{t) = exp {2\t) leads to the Gross-Pitaevskii equation of Ref. flG*! with 
an expulsive harmonic potential. Simulating the time dependencies of the harmonic potential and the interatomic 
interaction of Eq. ^ may be difficult experimentally, but if we choose the parameters such that ci S> ai 3> 0^2, 
the interatomic interaction can be considered practically as constant, and the strength of the harmonic potential will 
be oscillating as sin (uit -\- 5). It should be noted, however, that even in this limit, the contribution of the positive 
cos{Lot + 5)'^ term will result in that the trapping potential will spend more time being an inverted parabola than being 
an upright parabola. Over a large enough time interval, this will lead to expelling the solitons away from the center 
of the trap. This is indeed the picture that we get when we plot the density of a single-soliton solution as shown in 
Fig.[TJ It is clear from this figure that the soliton is being expelled out from the center of the condensate at a; = 0. The 
oscillation in the trajectory of the soliton is due to the oscillating trapping potential. The discontinuous appearing of 
the soliton's peak density is due to the interaction with the background. This trajectory can be also extracted from 
the general solution Eq. ([7]) by considering the term (rc^e^"''"^'^'-^ + Tc\e~'^°'^'^^^''^) in the denominator. At the 
peaks of the oscillation in Fig. ([1]), this is the dominant term that determines the position of the peak of the soliton. 



.dilj{x, t) 



dx'^ 



sin [iot + 5) + —a2 cos {ivt + S)^ 
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Specifically, the peaks are given by the condition 2ar] — v^A,.a; = 0. Using this condition to plot x vs. t in Fig. [2l 
we obtain a curve that is identical to the soliton trajectory in Fig. [TJ The mean slope of this curve is proportional to 
a/2Ar. Hence, the rate at which the soliton leaves the center of the trapping potential can be delayed by choosing 
the parameters and the arbitrary constants such that a/2Ar is small. 

For the special case of a = the center of mass of the soliton will be located at x = at all times. The oscillating 
trapping potential will result only in oscillations in the width and peak density of the solitons. The solitons in this case 
remain trapped at the center of the trapping potential. We show this case in Figs. [3][6] where we see a multi-soliton 
solution with central soliton being pinned at a; = and off-central ones oscillating around their initial positions. 
In Fig. m we use the same parameters as in Fig. |3] but with a doubled frequency lo. This shows clearly that the 
oscillations of the off-centered solitons are due to the oscillation in the trapping potential. In Fig. [5l we show that for 
some values of the parameters the dynamics of the peak soliton density can be drastic such that the soliton disappears 
in the background and reappears at regular discrete times. Figure [5] is similar to Fig. [3] but with a larger number of 
solitons. 



V. CONCLUSIONS 

We have found exact single and multi solitonic solutions of a time-dependent Gross-Pitaevskii equation with time- 
dependent amplitudes of harmonic trapping potential and interatomic interaction. We considered the interesting 
special case of sinusoidally oscillating strength of the harmonic trapping potential. We focused on the effect of the 
frequency and amplitude of the oscillating trapping potential on the dynamics of the solitons. The parameters can 
be chosen such that the soliton is trapped at the center of the condensate. For a typical ^^Rb condensate with 10^ 
atoms and trapping frequency of order 100 Hz, the Thomas-Fermi size of the condensate is Rtf = {Na/SY^^S w lOS 
[23 |. This means that the width of the solitons of Figs. [5][5]are of the order of the size of a nonsolitonic condensate 
while the width of the solitons in Figs[5][6]is smaller than the size of the nonsolitonic condensate. 
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FIG. 1: (Left) Spatiotemporal contour plot of solitons density profile. (Right) Soliton density profile aX t = 0. The values of 
the parameters used in this plot are: — C2 = ci = 10, 03 = 04 = 1, Aii = 2, Air- = 1, A = 0.1, a = 0.2, fco = 0, A = 1, w = 0.1, 
<J = 0, Qi = -4, aa = 0.3. 
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FIG. 2: Trajectory of the soliton peak density. The values of the parameters used here are the same as those of Fig.l. 
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FIG. 3: (Left) Spatiotemporal contour plot of solitons density profile. (Right) Soliton density profile at t = 0. The values of 
the parameters used in this plot are: — C2 = ci = 10, cs = C4 = 1, Aij = Air = 1, A = 2, a = 0.9, ko ~ 5.03, \ = uj = 1, 5 = 0, 
ai = —6, a2 = 0.3. The value of ko is the solution of a = with respect to fco. 
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FIG. 4: Spatiotemporal contour plot of density solitons density profile. The values of the parameters used in this plot are the 
same as those of Fig. 3 but with oj = 2. 
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FIG. 5: (Left) Spatiotemporal contour plot of solitons density profile. (Right) Soliton density profile at t = 0. The values of 
the parameters used in this plot are: — C2 = ci = 10, 03 = 04 = 1, Ai; = 0, Ai,. = 1, A = 2, a = 0.9, fco = 0, A = 1, a; = 0.01, 
5 = 0, ai = -2, 02 = 0.3. 
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FIG. 6: (Left) Spatiotemporal contour plot of solitons density profile. (Right) Soliton density profile at t = 0. The values of 
the parameters used in this plot are the same as those of Fig.3 but with ai = —1. 



